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We show that initially localized and uncorrelated two-particles quantum wavepackets evolving in 
a one-dimensional discrete lattice become strongly entangled while drifting under the action of an 
harmonic AC field resonant with doubled Bloch oscillations promoted by a static DC field. Although 
partial entanglement is achieved when the AC field is resonant with the single-particle Bloch oscil¬ 
lations, it is strongly limited by the survival of anti-correlated unbounded states. We further show 
that the phase dependence of the wavepacket centroid velocity is similar to the semiclassical behav¬ 
ior depicted by a single-particle. However, the drift velocity exhibits a non-trivial non-monotonic 
dependence on the interaction strength, vanishing in the limit of uncorrelated particles, that unveils 
its competing influence on unbounded and bounded states. 

PACS numbers: 03.67.Bg, 03.75.Lm, 37.10.Jk,, 67.85.-d 


Recent experiments investigating the behavior of few 
interacting atoms in optical lattices [l| opened the possi¬ 
bility to generate a variety of quantum entangled mat¬ 
ter states that can have a great impact on the search 
for universal and more efficient quantum computation 
processes, quantum sensing and metrology [1-Q. It has 
been experimentally demonstrated that quantum entan¬ 
gled atom pairs perform Bloch oscillations (BOs) at 
twice the fundamental frequency in tilted optical lat¬ 
tices (TOLs)[Tj. BO is the phenomenon of oscillatory 
motion of wavepackets placed in a periodic potential 
when driven by a constant forced, 0. It has been ob¬ 
served in several physical contexts such as semiconduc¬ 
tors superlattices 3, p , ultra-cold atoms [ol-fl^. optical^ 
and acoustic waves |14|. The coherent phenomenon of 
frequency doubling of BOs was predicted to occur in the 
presence of interaction (l5l - [l^ because two particles bind 
together and perform correlated tunneling. Fractional 
BOs at multiples of the fundamental BO frequency have 
also been demonstrated for larger clusters of interact¬ 
ing particles [l^. A photonic realization of the BOs fre¬ 
quency doubling has been proposed [2l|, and experi¬ 
mentally achieved in waveguide lattices (23|. Its experi¬ 
mental demonstration with ultra-cold atoms of bosonic 
®^Rb in decoupled one-dimensional tubes of an optical 
latticed represents an important step towards the in¬ 
vestigation of essential features of quantum many-body 
states. 

Unidirectional transport and super-BOs can be in¬ 
duced when wavepackets are driven by superposed static 
and harmonic fields [11-113. It is achieved when the har¬ 
monic field is resonant with the underlying BO frequency, 
which depends linearly on the strength of the static field. 
A small detuning from the resonant condition results 
in super-BOs due to an effective tunneling renormal¬ 
ization. This phenomenology has been experimentally 
demonstrated in a weakly interacting Bose-Einstein con¬ 
densate of Cs atoms placed in a TOL under forced driv¬ 
ing, achieving matter-wave transport over macroscopic 
distances [ 13 . It has been theoretically demonstrated 


that some features of super-BOs and unidirectional trans¬ 
port require an important phase correction to be include 
in the tunneling renormalization picture[28j|. Further, it 
has been theoretically shown that correlated super-BO of 
a bounded two-particles state can be observed under ap¬ 
propriate drive conditions [13 • Considering the current 
stage of experiments are in a position to probe the effect 
of interactions on driven transport of few correlated cold 
atoms [ 13 , it is of fundamental importance to have new 
schemes devised to generate and manipulate atomic en¬ 
tangled states in discrete lattices which can be used to 
probe quantum aspects of mater waves (lll - fsHl - 

In this letter, we show that two-interacting particles 
placed in a TOL can be coherently transported when an 
external harmonic field is made resonant with the dou¬ 
bled BO frequency. We also unveil the dependence of 
the wavepacket centroid velocity on the strength of the 
interaction as well on the relative phase of the harmonic 
field. Contrasting with the corresponding unidirectional 
transport for the resonant condition at the fundamental 
BO frequency, the degree of entanglement will be shown 
to continuously increase, with the particles developing 
positive spacial correlations due to the suppression of un¬ 
bounded wavepacket components in the entire range of 
reachable drift velocities. 

The dynamics of two interacting particles placed in 
a linear discrete lattice of spacing d under the action 
of superposed static DC and harmonic AC fields can be 
described in the framework of the tight-binding Hubbard 
model Hamiltonian as 

H— 'y ^ -|- + sF{t)dnb\^^bn^cr 
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where bn,a and 5l g. are the annihilation and creation 
operators for the particles of charge e at site n in spin 
state cr, J is the hopping amplitude, and U is the on-site 
Hubbard interaction. We are considering the particles 
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distinguishable by their spin state. In what follows we 
will use units ote = d= J = ti — 1. The applied field 
is assumed to be given by F{t) = Fq + cos (wt + (j)) 
where Fq stands for the amplitude DC field and F^ the 
amplitude of the AC field with frequency to. In the ab¬ 
sence of interaction, unidirectional transport occurs when 
the harmonic field is resonant with the fundamental BO 
frequency ujb = Fq and its sub-multiples . 

In order to follow the time evolution of the 
two-particles wavepacket, we solved the time de¬ 
pendent Schrodinger equation by expanding the 
wavevector in the Wannier representation |'I'(t)) = 

Eni.na V'ni.n 2 (i)l«- 1 :»^ 2 ), where the ket \ni,n 2 ) repre¬ 
sents a state with the particle in the spin state cr = I 
at site Til and the other particle at site 712 - The coupled 
motion equations for the wavepacket amplitudes 


= V’m + l.nsW +' 0 ni-l,77,2(0 +V'ni.n2-nW 

+ 0ni,n2-l(0 + F{t)ini + n 2 ) + ^ni ,712 C/Uni ,n2 (0 (2) 


were solved numerically by using a high-order method 
based on the Taylor expansion of the evolution operator 
A{6t) = exp {—iH5t) = , where FI is the 

two-particles Hamiltonian. Our calculations were taken 
by using St = 10“^ and the sum of evolution operator 
truncated at Uo = 12. The chain size was large enough 
to keep the wavefunction amplitude at the borders always 
smaller than I0“^°. We followed the time-evolution of an 
initial wavepacket composed of a non-entangled Gaussian 
wavepacket 


{ni,n2\'f{t = 0 )) = le-[Oi-"?)'+("2-n^)"]/4<7= ( 3 ) 

centered at the initial position (rii^n^)- A is a normal¬ 
ization factor. In the following numerical results, we will 
consider initial wavepackets with E = = 10 and 20 

and held strengths Fq = 0.6 and = 0.8Fo[2^, [l^. It 
is important to stress that no net displacement of the 
wavepacket centroid is achieved for an initial Fock state 
with the particles occupying a single siteQ. In this case, 
the wavepacket dynamics is symmetric with respect to 
the initial position. 

We start following the time evolution of the wavepacket 
centroid associated with each particle dehned as 

(n*(t)) = ^ (ni)|V'„i,„2(0P , i = l,2. (4) 

711,712 


Due to the symmetry of the initial state and interaction 
Hamiltonian, one has that (ni(t)) = (n 2 (t)). In hg. 1 
we plot the centroid evolution for both cases of the AC 
held resonant with the fundamental and doubled BO fre¬ 
quencies for an intermediate interaction strength. A net 
unidirectional modulated motion of the wavepacket cen¬ 
troid is obtained, whose average velocity depends on the 
relative phase of the AC held. 



FIG. 1: Time evolution of the one-particle wavepacket cen¬ 
troid. The initial state is composed by a direct product of 
Gaussian wavepackets with E = 20 centered at {ni,n 2 ) = 
(0,0). The interaction strength is H = 3. The AC and DC 
fields satisfy the resonance conditions u = Fq (one-particle 
BO frequency) and lj = 2Fo (correlated BO frequency). Re¬ 
sults for distinct phases <j) of the AC held are shown to stress 
the phase-dependence of the asymptotic centroid velocity. 


As shown in hg. 1, a phase control of the AC held allows 
to tune both the speed and direction of the centroid mo¬ 
tion. Within a semiclassical description, the wavepacket 
net velocity of a non-interacting particle driven by an AC 
held in resonance with the BO is given by [ 2 ^ 

V = 2Ji{Fuj/Fo)cos[{Fuj/Fo)cos(j) - </>] (5) 

where Ji{x) is the Bessel function of the hrst kind of 
order 1. In hg. 2, we plot the phase dependence of 
the average centroid velocity for the case of interacting 
particles driven by a resonant AC held. The phase de¬ 
pendence converges to the above semiclassical prediction 
as the initial wavepacket becomes wider, although pre¬ 
senting distinct overall amplitudes and reversed trends 
for the cases of resonance with the fundamental and 
doubled BO frequencies: larger (smaller) velocities are 
reached for wider initial wavepackets at the fundamen¬ 
tal (doubled) BO resonance. No net transport occurs 
for (j) = ±7r/2 while maximum speeds are reached at 
4>i = {Fuj/F q) cos <j)i ~ 0.64II... (negative velocity) and 
(j )2 ='^ — 4’! = 2.5008... (positive velocity). 

Although the phase dependence of the drift velocity is 
similar to the one depicted by non-interacting particles 
driven by an AC held resonant with the usual BO fre¬ 
quency, the drift velocity at the doubled BO frequency 
is, ultimately, an interaction effect. In the absence of 
interaction {U = 0), there is no net transport at the dou¬ 
bled BO resonance. Unidirectional transport of a single 
particle can only be promoted when the frequency of the 
AC held is the fundamental BO frequency (or its sub¬ 
multiples) [H, [ 2 ^. 

Considering that the actual value of the inter-particle 
interaction can be tuned via the depth of optical 
lattices [I|, we computed the extremal values of the cen¬ 
troid velocity as a function of U for both resonant cases 
(see hg. 3). The numerical result at C/ = 0 and uj = Fq 
is consistent with the semiclassical prediction Vmax = 
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FIG. 2: Phase dependence of the one-particle centroid veloc¬ 
ity for U — 3, two distinct initial wavepackets, uj = Fq (left 
panel) and lj = 2Fo (right panel). The solid lines correspond 
to the semiclassical dependence v oc cos (-Ss^ cos((()) — ((i) to 
which the results converge as E —>■ oo. Notice the reversed 
dependence of the maximum velocity with E at the single 
particle and correlated particles resonances. 


2Ji(J1j/^o = 0.8) = 0.73768... for non-interacting par¬ 
ticles. When the interaction is turned on, the drift cen¬ 
troid velocity at the fundamental resonance initially de¬ 
creases, passing at a minimum value on an intermediate 
coupling strength. On the other hand, a net transport 
develops at the doubled BO resonance, reaching a maxi¬ 
mum also at a finite U. The maximum velocity reached 
at the correlated two-particles doubled BO resonance is 
of the same magnitude of the minimum drift velocity at 
the fundamental BO resonance. Such non-monotonous 
dependence of the drift centroid velocity on the cou¬ 
pling strength is due to two opposite effects played by 
the interaction in the wavepacket dynamics. The initial 
wavepacket can be viewed as composed of components 
associated with bounded and unbounded states. The un¬ 
bounded wavepacket components are not drifting under 
the action of the frequency doubled AC field. Therefore, 
the observed unidirectional motion of the wavepacket 
centroid is solely due to the drifting of the bounded 
wavepacket components. In these bounded states, the 
system behaves as a composite particle of charge 2s. 
Therefore, while the velocity dependence on the field fre¬ 
quency is the same as that of a single (although com¬ 
posite) particle, its non-monotonic dependence on the 
interaction strength U unveils its distinct influence on 
bounded and unbounded states. For the components as¬ 
sociated with bounded two-particle states, the interac¬ 
tion favors the pairing and promotes coherent hoppings. 
For the unbounded components, the interaction enhances 
the wavepacket width, thus reducing the double occu¬ 
pancy probability. As a result of these competing ef¬ 
fects, optimal correlated two-particles motion occurs at 
an intermediate coupling strength[l^. It has been previ¬ 
ously reported that such competition also leads to a non¬ 
monotonic behavior of the BO amplitudefl^ and Ander¬ 
son localizationfH, 113 . 

To get more physical insight on the wavepacket dy¬ 
namics at the fundamental and doubled BO resonances. 



FIG. 3: Maximum one-particle drift velocity as a function 
of the interaction strength U for an initial wavepacket with 
E = 20, ui = Fo (left panel) and ui = 2fo (right panel). The 
solid lines are guides to the eyes. The interaction suppresses 
the maximum drift velocity when the AG field is resonant with 
the single-particle BOs, reaching a minimum at an intermedi¬ 
ate value of U. The counterpart peak in the maximum drift 
velocity for the resonant at the doubled BO frequency reflects 
the predominant role played by the two-particles bounded 
states in this regime. 


we plot in fig. 4 some two-particles wavepackets after 
a finite evolution time. Three representative phases of 
the AC held were chosen, corresponding to zero, inter¬ 
mediate and maximal drift velocities. The interaction 
strength {7 = 3 favors the coherent two-particles mo¬ 
tion. At the fundamental BO resonance, the wavepacket 
develops well distinct structures. The one along the di¬ 
agonal stands for bounded states with the particles lo¬ 
cated spatially close. The structure away from the di¬ 
agonal stands for unbounded states with the particles 
driven to opposite sides of the chain. Only near the phase 
leading to maximal drift velocity the unbounded compo¬ 
nent is suppressed. When the AC held is resonant with 
the correlated two-particles doubled BO frequency, the 
wavepacket spreads while keeping the particles spatially 
close for all values of the AC held phase. Under this con¬ 
dition, the wavepacket is composed mainly of bounded 
two-particles states. This feature indicates that strong 
quantum correlations are always developed between the 
two-particles over several lattice sites. Such spatially ex¬ 
tended quantum entangled state can be explored to probe 
essential quantum aspects of matter waves, such as non¬ 
locality. 

There are several prescriptions to quantify the degree 
of entanglement of a two-particles wavefunction. Here, 
we will use as a diagnostic tool the purity function de- 
hned as F(t) = Trp\(t), where pi(t) is the reduced den¬ 
sity matrix for particle 1 obtained after taking the partial 
trace over the states of particle 2 {pi{t) = Tr 2 p{t), with 
p{t) = |4'(t))(4'(t)|). From fundamental properties of the 
density matrix, the purity function P{t) = 1 for a pure 
state, meaning that the two particles are not quantum en¬ 
tangled. It assumes the value P{t) = 1/N whenever the 
quantum state of particle 1 is an even incoherent distribu¬ 
tion among N states. In the present scenario, when the 


















4 









200 



200 




m 



too 





% 


0 


/ 


- too 



-too 




-200 

10= 1.57 


-200 

= 

i . r )7 



-200 -100 0 100 200 



-2<U> 

- too 0 too 

200 








200 



200 




100 



too 





V 1 

^1 

0 


/ 



. } 






-im 



- too 




-200 

0 ^ 2.00 


-200 

<0 = 

2.00 



-200 - too 0 too 200 



-200 

-m) 0 too 

200 








200 



200 




m 



too 





/ 


0 




-m 



- too 




-200 

0 = 2 ..50 


-2tMt 

<t> — 

2.50 



-200 - too 0 loo 200 



-200 

-tm 0 m 

200 









FIG. 4: (Color online) Density plots of wavepacket after a 
finite evolution time and interaction strength [7 = 3. Top 
panels: (j) = 7r/2 for which the centroid velocity is zero; In¬ 
termediate panels: </> = 2.0 leading to an intermediate ve¬ 
locity; Bottom panels: 0 = 2.5 resulting in maximal veloc¬ 
ity. Left panels: AC field resonant with the single-particle 
BO (cj = Fq). The wavepacket develops positively corre¬ 
lated components (structure along the diagonal representing 
bounded particles) as well as anti-correlated components with 
the particles being driven to opposite sides (unbounded par¬ 
ticles), except near the maximal drift velocity. Right panels: 
AC field resonant with the correlated-particles BO (a; = 2Fo). 
The wavepacket remains along the diagonal (bounded parti¬ 
cles) irrespective to the phase. 


two particles become maximally quantum entangled over 
N lattice sites, i.e., Idf) = {1/\/n) |1 : |2 : m), 

the partial density matrix pi = {^/N)J2^=i |1 • ^i)(l ^ 
nil the purity function reaches P{t) = 1/N. In terms 
of the wavefunction components, the purity function can 
be written as 

P ~ 'y ( 6 ) 

ni ,mi ,712 l'^2 

As a complimentary tool, we also computed the two- 
particles normalized correlation function defined as 

C{t) = [(nins) - (ni)(n2)]/[(ni)(n2)]- (7) 
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FIG. 5: Time evolution of the quantum purity function (left 
panels) and two-particles correlation function (right panels) 
for distinct phases of the AC held and U = 3. (a-b): AC 
held with frequency ui = Fq. (c-d): AC held with frequency 
uj = 2Fo. At w = 2Fb, the predominance of bounded states 
generates positive correlations and leads to an increased de¬ 
gree of quantum entanglement (continuously decreasing pu¬ 
rity measure) irrespective to the phase of the AC held, (e-f): 
The cases of initial Fock states (E = 0) with particles occupy¬ 
ing the same do = 0 and well separated do = 10 sites and held 
phase (f) = 2.5. Although there is no net transport, a continu¬ 
ously increasing degree of quantum entanglement and positive 
correlations is only obtained when initially close particles are 
driven at to = 2 Fq. 


In fig.5 we show the time evolution of the purity and 
pair-correlation functions in the presence of and AC field 
resonant with the fundamental and doubled BO frequen¬ 
cies. When the particles are driven by an AC field at 
resonance with the fundamental BO frequency, the de¬ 
gree of entanglement saturates (although slower as the 
extremal drift velocities are approached). Such satura¬ 
tion is related to the fact that unbounded states pre¬ 
ponderate in the wavepacket dynamics (see fig. 4). The 
component with particles driven to opposite sides of the 
chain has strong anti-correlations which can surpass the 
positive correlation associated with the bounded compo¬ 
nents. Therefore, the net pair-correlation results nega¬ 
tive in a finite range of phase values. In contrast, the 
purity function continuously decrease in time when the 
particles are driven at resonance with the doubled BO 
frequency, irrespective to the phase of the AC field. This 
feature indicates that the wavepacket develops quantum 
entanglement over a continuously growing chain segment 
for the entire range of allowed drift velocities. The pos¬ 
itivity of the pair-correlations reflects the predominant 
role played by bounded two-particles states. In Fig.5(e- 
f) we report results for the case of initial Fock states 
with the particles occupying the same site, as well as 
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separated sites. Even depicting no net displacement of 
the wavepacket centroid, continuously increasing quan¬ 
tum entanglement and positive pair-correlations are still 
obtained for w = 2 Fq when the particles initially occupy 
the same site. 

To conclude, we recall that recent experiments on cold 
atoms trapped in TOLs showed the capability of observ¬ 
ing coherent atomic motion under the action of an AC 
field over macroscopic distancesfl^ Further, the co¬ 
herent BO of two quantum entangled particles has also 
been demonstrated experimentally [l|. Therefore, the 
presently proposed scheme to generate and manipulate 
spatially extended entangled two-atoms states by driven 
them using an AC field resonant with the frequency dou¬ 
bled two-particles BO is well within currently accessible 
experimental techniques. The speed and direction of the 
wavepacket centroid motion can be externally controlled 
by the AC field while the strength of the interaction can 
be changed by tuning the potential depth. While the 
phase dependence of the drift velocity is similar to the 
one displayed by non-interacting particles driven by an 


AC field resonant with the usual BO frequency, its de¬ 
pendence on the interaction strength U is rather non¬ 
trivial. The drift velocity vanishes in the limit of non¬ 
interacting particles and varies non-monotonically with 
U, unveiling the distinct roles played by unbounded and 
bounded states components on the wavepacket dynam¬ 
ics. Very recently, entanglement measures have been used 
to experimentally characterize the dynamics of strongly- 
correlated many-body systems [s^. Although decoher¬ 
ence effects due to inherent coupling with other degrees 
of freedom shall be carefully taken into account, the real¬ 
ization of spatially extended entangled two-atoms states 
by exploring the phenomenon of frequency doubling of 
BOs may impel the development of a new class of ex¬ 
periments aiming to search for signatures of quantum 
non-locality in matter waves. 
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